Deep Mixing in Evolved Stars: 
I. The Effect of Reaction Rate Revisions from C to Al. 
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ABSTRACT 

We present computations of nucleosynthesis in low-mass red-giant-branch and 
asymptotic-giant-branch stars of Population 1 experiencing extended mixing. We 
adopt the updated version of the FRANEC evolutionary model, a new post- 
process code for non-convective mixing and the most recent revisions for solar 
abundances. In this framework, we discuss the effects of recent improvements 
in relevant reaction rates for proton captures on intermediate-mass nuclei (from 
carbon to aluminum). For each nucleus we briefly discuss the new choices and 
their motivations. The calculations are then performed on the basis of a pa- 
rameterized circulation, where the effects of the new nuclear inputs are best 
compared to previous works. We find that the new rates (and notably the one 
for the ^^N(p,7)^^0 reaction) imply considerable modifications in the composi- 
tion of post-main sequence stars. In particular, the slight temperature changes 
due to the reduced efficiency of proton captures on ^^N induce abundance varia- 
tions at the first dredge up (especially for ^^O, whose equilibrium ratio to ^^O is 
very sensitive to the temperature). In this new scenario presolar oxide grains of 
AGB origin turn out to be produced almost exclusively by very-low mass stars 
(M < 1.5 — 1.7 Mq), never becoming C-rich. The whole population of grains 
with ^^O/^^O below 0.0015 (the limit permitted by first dredge up) is now ex- 
plained. Also, there is now no forbidden area for very low values of ^'^O/^^O 
(below 0.0005), contrary to previous findings. A rather shallow type of transport 
seems to be sufficient for the CNO changes in RGB stages. Both thermohaline 
diffusion and magnetic-buoyancy-induced mixing might provide a suitable physi- 
cal mechanism for this. Thermohaline mixing is in any case certainly inadequate 
to account for the production of ^^Al on the AGB. Other transport mechanisms 
must therefore be at play. In general, observational constraints from RGB and 
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AGB stars, as well as from presolar grains, are well reproduced by our approach. 
An exception remains the nitrogen isotopic ratio in mainstream SiC grains. For 
the low values measured in them (i.e. for ^''N/^^N < 2000) we have no expla- 
nation. Actually, for the several grains with subsolar nitrogen isotopic ratios no 
known stellar process acting in low mass stars can provide a clue. This might be 
an evidence that some form of contamination from cosmic ray spallation occurs 
in the interstellar medium, adding fresh ^^N to the grains. 

1. Introduction 

It is now well known that the processes of astration and gas return to the Interstellar 
Medium (ISM) of galaxies is dominated by low mass (LM) and intermediate mass (IM) 
stars, i.e. by objects below M = 7-8 Mq (Sedlmayr 1994). These stars share the properties 
of ascending, after the Main Sequence, along the Red Giant Branch (RGB) and then, after 
the occurrence of core-He burning, of terminating their evolution through the so-called 
Asymptotic Giant Branch (AGB) stage (sec e.g. Iben & Renzini 1983). Here they lose mass 
efficiently thanks to radial pulsations and to radiation pressure on dust grains (Herwig 
2005; Guandalini et al. 2006; Guandahni & Busso 2008), powering stellar winds that are 
enriched with the products of the internal nucleosynthesis. Tiny solid particles that are 
abundant in such winds maintain the signature of those nuclear processes; samples of the 
most refractory among those grains can be recovered from pristine meteorites (Zinner 2005) 
and this fact has enormously increased the importance of AGB stars as laboratories for the 
study of H- and He-burning nucleosynthesis (Busso et al. 1999). 

At the end of the AGB stage, the residual envelopes of LM and IM stars are lost 
through the fast ejection of hot, ionized materials ("superwind"); the circumstellar envelope, 
heated and illuminated by the now hot central star, forms a planetary nebula, while the 
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dense, electron-degenerate core starts a blue-ward path, which will ultimately give birth 
to a white dwarf (Van Winckel 2003; Prada Moroni et al. 2004). Stars crossing these 
evolutionary stages not only represent a fundamental component of the baryonic mass in 
galaxies, but also provide us with invaluable information on the physical and chemical 
evolution of the universe. As a consequence, the understanding of moderately massive stars 
is a necessary condition for addressing a large part of modern astrophysical problems. 

Unfortunately, embarrassing gaps exist in our knowledge of their physics, and 
fundamental problems affecting their evolution are not properly treated in common stellar 
modeling. This is not surprising, as stellar models provide only a basic description, 
where important physical mechanisms (like convection) are oversimplified and others (e.g. 
rotation) are often neglected (see e.g. Charbonnel & Lagarde 2010, for a recent outline of 
the importance of rotation). One big problem affecting low- mass-star evolutionary models 
is revealed by the isotopic ratios of CNO elements, as observed in stellar photospheres 
(Harris et al. 1985; Cottrell & Sneden 1986; Gilroy & Brown 1991; Shetrone et al. 1993; 
Smiljanic et al. 2009) and in presolar grains of AGB origin, preserved in meteorites (see 
Nittler et al. 1997; Choi et al. 1998; Amari et al. 2001). 

Some chemical peculiarities concerning light elements are found very early in stellar 
evolution (Richard et al. 2005; Melendez et al. 2009). Others appear shghtly after the 
occurrence of the first dredge up (hereafter FDU) , where stars populate the so-called bump 
of the luminosity function (Charbonnel 2004). Then the appearance of anomalies in the 
composition continues along the AGB phase (Wasserburg et al. 1994; Busso et al. 2010). 
The observational basis for these issues has grown considerably over the years, confirming 
chemical anomalies at least from Li to O, but sometimes extending up to Mg and Al, 
including the production of the unstable isotope ^®A1 (sec e.g. Pilachowsky et al. 1993; 
Charbonnel 1994; Gratton et al. 2000; Grundahl et al. 2002; Busso et al. 2003; Zinner et al. 
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2006; Nittler et al. 2008). 

The above evidence suggests that non-convective transport mechanisms occur in 
stars, especially during the RGB and AGB phases, linking the envelope to regions where 
proton captures occur (Sweigart & Mengel 1979; Boothroyd et al. 1994; Charbonnel 1995; 
Charbonnel & Do Nascimiento 1998; NoUett et al. 2003) 

Secular mixing phenomena are made easy in RGB and AGB phases immediately 
following an episode of convective dredge-up, because the H-burning shell then advances in 
homogeneous regions, so that the natural barrier opposed to mass transport by the chemical 
stratification is not present. Mechanisms of this kind might actually account for several 
puzzling astrophysical observations: from the evolution of ^He in the ISM (Wasserburg et 
al. 1994; Glocker & Geiss 1998; Bania et al. 2002) to the complex phenomena of production 
and destruction of Li in the Galaxy, including its long-term decline in main sequence 
stars (Spiegel 1972; Spiegel and Zahn 1992; Zahn 1994) and its re-production in some 
(rare) red giants (Charbonnel & Balachandran 2000; Guandalini et al. 2009). Further 
implications of slow mixing probably include the anti-correlations displayed by globular 
cluster stars, involving intermediate mass elements up to Mg and appearing already on the 
Main Sequence (Gratton et al. 2001; Thevenin ct al 2001; Ramirez & Cohen 2002). Indeed, 
they might have been inherited from LM or IM progenitors, hosting proton captures in 
mixing episodes through the AGB phases. 

Among the processes involved in deep-mixing phenomena, those induced by rotation 
certainly play an important role (Zahn 1992; Charbonnel & Lagarde 2010), despite 
difficulties in explaining the observed isotopic ratios (Palacios et al. 2006). Actually, a 
shear layer, at the contact between the differentially rotating convective envelope and the 
more rigid stellar core is a natural expectation (see e.g. Kippenhahn & Weigert 1990). 
A rigid-body behavior in the cores of low-mass stars might be linked to an interplay of 
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magnetic coupling and meridional circulation (Eggenberger et al. 2005). 

Despite formal differences in the various approaches, most models considered the 
chemical mixing induced by rotational effects (and the associated angular momentum 
transport) as a diffusive process, with parameterized diffusion coefficients (see e.g. 
Denissenkov et al. 1998; Denissenkov & Weiss 1996, 2000; Weiss et al. 2000). This approach 
has also been adopted in modeling massive, stratified stars, e.g. by Maeder & Meynet 
(2004). Alternatives were presented by Denissenkov & Tout (2003), who studied the effects 
of gravitational waves, and by Busso et al. (2007a); Nordhaus et al. (2008); Denissenkov 
et al. (2009). These last authors discussed the effects of magnetic buoyancy induced by a 
stellar dynamo. The presence of such a dynamo in red giants finds supports in some stellar 
variability observations (Andrews et al. 1988) and in the identification of surface magnetic 
fields in AGB stars (Herpin et al. 2006). Stellar magnetic fields and their instabilities have 
been also the object of previous detailed analysis (Spruit 1999, 2002). However, by far the 
most popular scenario in recent years has been that of the so-called thermohaline diffusion, 
induced by an inversion in the molecular weight yU. In H-burning regions this is due to the 
activation of the reaction ^He+^He ^He + 2p (Eggleton et al. 2006, 2008; Charbonnel 
& Zahn 2007), which locally reduces The phenomenon is more general than that; it 
was recognized long time ago and it can occur in other astrophysical contexts (Stothers & 
Simon 1969; Ulrich 1972). 

Diffusive mechanisms (such as thermohaline mixing) and circulation phenomena (like 
those induced e.g. by magnetic buoyancy) yield abundance changes that are sometimes 
undistinguishable from one another, sometimes remarkably different, depending on the 
velocity of the transport. Verifying whether observations can provide clues to discriminate 
among the possible physical processes generating the mixing is therefore an important 
issue. In this respect, interesting new elements have been offered by the recent 2D- and 
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3D-calculations by Denissenkov (2010); Denissenkov & Merryfield (2010). After presenting 
our results, we shall comment in some detail on these elements, especially in comparison 
with what magnetic buoyancy can a priori provide. 

Before the effects of different, uncertain mixing models can be addressed, however, 
one has preliminarily to ascertain that the normalizing solar abundances and the nuclear 
input data be very reliable. Using all the upgrades appeared in the recent literature on 
these issues, we intend to revisit the coupled processes of proton capture nucleosynthesis 
and of deep mixing in red giants, starting from a parameterized analysis (where nuclear 
and abundance effects are best compared to previous work). In section 2 we shall present a 
general discussion of the set of new inputs adopted, postponing to a dedicated Appendix the 
details of the reaction rate choices made. In section 3 we shall then discuss recent upgrades 
and novelties in stellar models and in section 4 we shall briefly illustrate our technique 
for performing parameterized, post-process deep-mixing calculations, starting from the 
results of the FRANEC evolutionary code. Then in sections 5 and 6 we shall outline 
the nucleosynthesis results thus obtained, showing some comparisons with previous work, 
with stellar observations and with presolar grain measurements of isotopic abundances for 
elements from C to Al. Comments on recent models for thermohaline mixing in 2 and 
3 dimensions are then added in a dedicated section (section 7). Finally, some general 
conclusions are drawn in section 8. 

2. On Recent Revisions in Soleir Abundances and Reaction Rates 

A previous parameterized analysis of extra-mixing in red giants (NoUett et al. 2003) 
adopted solar abundances from Anders & Grevesse (1989), coupled with the general 
prescriptions for reaction rates from the compilation by Angulo et al. (1999), hereinafter 
simply referred to as NACRE. Since then, however, many remarkable changes have occurred 
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in both sets of input data. For what concerns solar abundances, we shall in general adopt 
the compilation by Asplund et al. (2009), implying a total solar metallicity of Z = 0.014. 
We underline that not only elemental, but also isotopic abundances (as updated in that 
review) are important for our analysis; this is especially the case for CNO elements. In 
particular, we shall make use of the new estimate of the ^^N/^^N ratio, as deduced by 
studies of the Jovian atmosphere (Fouchet et al. 2000; Owen et al. 2001); for it we adopt 
the value derived by the Composite Infrared Spectrometer onboard the Cassini spacecraft 
(2.22±0.52xl0~^, a factor 1.5 smaller than the terrestrial one), which should represent 
the initial solar-system ratio (Fouchet et al. 2004). For oxygen isotopes we use the values 
suggested by Asplund et al. (2009). 

Concerning reaction rates, recent extensive reviews have been published, summarizing 
the main experimental changes occurred after the NACRE compilation, which has been 
the reference in most papers on low-mass stellar nucleosynthesis over the last decade. In 
particular, we shall discuss here (and adopt throughout the paper) the recommendations by 
Adelberger et al. (2010) and lliadis et al. (2010), hereinafter simply refereed to as ADEIO 
and ILIIO, respectively. We instead will not analyze the work by Beard et al. (2010), 
which pertains to stellar burning at much higher temperatures and to nucleosynthesis in 
high-density environments. The work by ADEIO represents the follow-up of the solar-fusion 
cross-section review by Adelberger et al. (1998) and includes subsequent measurements for 
reactions of the pp chain and of CNO cycles. It mainly focuses on low-energy cross sections, 
aiming at extracting the rates for temperatures typical of the solar core (Tg^^ 0.015). 
Reaction rates are obtained analytically, in terms of low-cncrgy expansion coefficients of 
the astrophysical S(E) factors. The work by ILIIO is instead an upgrade of the previous 
review by lliadis et al. (2001). In the new work several improvements have been introduced: 



^For the temperature we adopt the notation where Tg =T/10^K. 
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(i) a wider range of atomic masses is covered, from A=14 to A=40, including some CNO 
reactions; (ii) updated nuclear physics inputs are considered in the reaction rate evaluation; 
(iii) a calculation procedure improved with respect to that of Iliadis et al. (2001) is 
adopted. In detail, a suited probability distribution is first associated with each nuclear 
physics quantity entering the reaction-rate calculation. Then, the reaction rate itself is 
computed for a set of input parameter values, which are randomly taken, each following the 
appropriate distribution. In this way, a probability density function is generated for the 
reaction rate reflecting the distribution of the input values. This method has the important 
consequence that the upper and lower limits are estimated with a statistically-appropriate 
approach, so that a really meaningful value can be attributed to them. Moreover, the rate is 
obtained over a wide temperature range Tg — 0.01 — 10, almost reaching down to the values 
found at the inner borders of the convective envelopes in low-mass AGB stars (Tg ~ 0.005). 
It is worth remarking that both compilations, ADEIO and ILllO, account for the recent 
results obtained using indirect techniques, such as the Coulomb Dissociation (Baur & Rebel 
1994), the Asymptotic Normahzation Coefficient (Mukhamedzhanov et al. 1997) and the 
Trojan Horse Method (La Cognata et al. 2010a). This makes the low-temperature reaction 
rates more accurate, thanks to the possibility of extending the measurements down to the 
astrophysical range of energies. 

The recommended reaction rates from ADEIO and ILIIO differ, sometimes remarkably, 
from NACRE results. This fact (and the search for the descending astrophysical 
consequences) is the main motivation for the present work. 

Proton capture reactions on nitrogen and oxygen isotopes are of special importance 
in our context. In particular, the reaction rate for ^^N(p,7)^^0 by ADEIO, reduced by 
50% with respect to NACRE, leads to a reduction in the efficiency of the whole CNO 
cycling, moving the H-burning shell toward inner stellar regions, where temperature and 
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density values are higher by about 10% and 25%, respectively, as compared to previous 
stellar models. We shall discuss at length the important consequences of these facts on the 
envelope abundances. As we shall see, also the new data for proton captures on ^^O induce 
non-neghgible effects on the surface composition. By contrast, no appreciable effect derives 
from the reduction in the rate for ^^N(p,7)^^0; in fact, it modifies the abundances of ^^O, 
^^O, ^^O and ^^F only in the regions where the H-burning efficiency is maximum (Tg > 
0.05), which are untouched by deep mixing. 

In general, we use the ADEIO compilation for pp-chain and CNO-cycle reactions, 
while we adopt the results from ILIIO for heavier nuclei. However, the two compilations 
partly overlap in the region of CNO isotopes. In this case, specific choices are necessary. In 
order to motivate them, and to illustrate some of the novelties with respect to NACRE, we 
present in the Appendix a discussion of the available data for some key reactions and of our 
choice. 

3. The new reference stellar models 

Although we have certainly achieved a general understanding of the conditions under 
which chemical elements are produced in stellar interiors, a predictive nucleosynthesis 
theory based on first principles is still lacking. A basic missing ingredient in developing 
such a theory is a full understanding of the transport properties of the turbulent, reactive 
plasma in which the elements are forged. Within this context, many efforts have been 
spent in the past years to study the properties of the final phases of stellar evolution by 
means of theoretical models calculated with the Frascati RAphson-Newton Evolutionary 
Code (FRANEC: see Chieffi et al. 1998). As an example, the extra- mixing calculations 
by NoUett et al. (2003), which represent the starting point of the present work, were 
based on FRANEC stellar structures taken from Straniero et al. (2003). This last work 
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presented parameterized theoretical formulae permitting a simple derivation of the physical 
characteristics of AGB stars (such as, for example, the number of thermal pulses, the mass 
of the H-exhausted core and of the envelope, the stellar surface temperature and so on), by 
interpolating the outputs of previously-computed full AGB evolutionary models (Straniero 
et al. 1997). 

Calculations presented in this paper are instead based on new evolutionary stellar 
sequences. In particular, we have performed complete runs (from the pre-Main Sequence to 
the end of the AGB phase) for a 1.5 Mq star and a 2.0 Mq star, with a metallicity equal 
to the most recent estimate for the solar one (Cristallo et al. 2009a). The runs have been 
performed with an upgraded version of the FRANEC code. Gomputations extending up 
to the RGB tip have also been included for various other masses in the range 1-1.7 Mq. 
Our calculations do not yet include rotation, nor the effects of magnetic fields. However, 
as compared to older versions of the FRANEG stellar code (Straniero et al. 1997), the 
models here adopted represent a remarkable improvement, including several new or updated 
ingredients. 

Among the main differences with respect to the previous scenario we recall: i) the use 
of a full nuclear network, coupled with the physical evolution; ii) the adoption of an updated 
solar chemical distribution (in particular, the code was revised in the last years using the 
Lodders 2003, compilation).^ We remind that older versions of FRANEG used instead the 
abundances by Anders & Grevesse (1989), which are substantially different; iii) the use 
of new reaction rates, of which crucial is in particular the new value for the ^''N(p,7)^^0 
reaction (see the Appendix for details); in the FRANEG code this update was originally 

^In the post process we use the updated suggestions by Asplund et al. (2009): the small 
differences between this set of data and the one used in the stellar code does not introduce 
any consequence on the model physics. 
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implemented when the measurements by the LUNA experiment became available (Imbriani 
et al. 2005; Lemut et al. 2006); iv) the adoption of new low-temperature opacities and their 
modifications during the evolution, when the envelope chemical admixture is changed by 
dredge up and becomes carbon rich (Cristallo et al. 2007); v) the choice, in evolved stages, 
of a mass loss law different from the one previously used; this relation is similar to the 
one by Vassiliadis & Wood (1993), but adopts more recent data on the Period-Luminosity 
relation of AGB stars (see Straniero et al. 2006, for details); vi) a different treatment of 
the radiative/convective interfaces. In particular, instead of using a bare Schwarzschild 
criterion, an exponentially-decreasing profile of the convective velocities was applied at the 
inner border of the convective envelope (Cristallo et al. 2009a). 

In any case, for a discussion of the improvements outlined above as items from i) to 
vi) and of their motivations see especially Cristallo et al. (2009a) . For the scopes of this 
work, we want primarily to underline that the adoption of a full nuclear network (although 
it increases the computational time by a factor 15 to 40, depending on the case considered) 
has proven to be crucial. Indeed, computing a set of full models, with the evolution of all 
the abundances coupled to the structural and thermodynamic evolution, provides bona-fide 
abundance templates, on which to verify the outcomes of any post-process calculation 
done for specific purposes. These post-process studies remain necessary whenever many 
parameterized cases have to be considered (as is required, e.g., in the present paper). 
However, now their results (in terms of predicted abundances) can be validated carefully on 
a smaller set of complete FRANEC calculations. The inclusion of the extended network into 
the code was primarily motivated by other reasons. In fact, the usual approach of stellar 
models, to consider only the limited set of reactions controlling the energy generation, 
becomes completely inadequate whenever the physical properties of a structure are strongly 
related to the chemical composition changes. This is, e.g., the case of the proton ingestion 
by the convective shell driven by the first fully developed thermal pulse, which is found to 
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occur in very metal-poor AGB stars and during the post-AGB phases at larger metaUicities. 
(Cristallo et al. 2009b; Herwig et al. 2010). More generally, the code with the full network 
can now be used to compute the nucleosynthesis in the He-shell (including s- processing). In 
the context of the present work verifying the post-process calculations with a (limited) set 
of checks from real, full-network models was important for those isotopes whose abundances 
are strongly influenced by the TDU. Typical examples of this are nuclei like ^^Ne (produced 
in the thermal pulses), ^^Na (that derives from proton and neutron captures on these same 
^^Ne seeds) and ^^Al. This last nuchde is produced in the innermost H-burning zones, hence 
sinks into the He-rich buffer and is subsequently largely destroyed by neutron-captures at 
high temperature (T > 3 x lO^K) during the thermal pulses. The surviving part is then 
dredged up again to the H-rich layers and can account for ^^Al/^^Al ratios in the envelope 
of up to 5 — 6 X 10""^, which are then enhanced by the extra production of ^^Al in deep 
mixing episodes. 

As we shall discuss in next section, the new models imply important changes in the 
surface isotopic composition already on the RGB, at the occurrence of the FDU. In the 
final evolutionary stages, then, the effects are even more substantial, especially for what 
concerns the chemical enrichment of the envelope along the AGB. 

4. Computations of Extended Mixing and Nucleosynthesis 

Starting from the models described in section 3, and also for the sake of comparison 
with previous works, we have performed parametric calculations, using as free parameters 
the rate of mass transport (M ) and the maximum temperature sampled by the circulating 
material (Tp), as in NoUett et al. (2003). This last parameter is best expressed by the 
difference A = logTf/ — log Tp, where is the temperature at which the energy from 
the H-burning shell is maximum. We shall use this definition of A throughout this paper; 
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similarly, we shall often express M in units of 10~^ Mq/jic and we shall call this parameter 
Mg- We present here results mainly for solar- metaUicity models of 1.5 and 2.0 Mq, although 
some extra models have been computed for comparisons with observational data. The 
post-process calculations have been run in a way similar to what was done previously, using 
older releases of the same FRANEC stellar code (NoUett et al. 2003; Busso et al. 2010). 
The post-process program is new and its results have been cross-checked with the ones 
obtained with the older code (Palmerini et al. 2008). 

During RGB stages, we have introduced extra- mixing after the "bump" of the 
luminosity function (BLF), i.e. after the H- burning shell has erased the chemical 
discontinuity left behind by the FDU (Charbonnel 1996). We shall not repeat here a 
general discussion of RGB physics, as it can be easily found elsewhere (Charbonnel & 
Balachandran 2000; Palacios et al. 2006). Our procedure in the calculations was illustrated 
by Palmerini et al. (2009). We have followed the abundance changes during the transport 
due to nucleosynthesis (partial derivatives with respect to time) and to displacements across 
the integration grid (partial derivatives with respect to distance, times the mixing velocity). 
One can write this synthetically as: 

dYi dYi ■ dYi 

— — - + M 

dt dt dM 

Here Y^ is the abundance by mole of the i-th isotope (for the present work, this index 
spans from H to ^^S). In agreement with NoUett et al. (2003) and with Busso et al. (2010) 
the fractional areas and fd occupied by the upward and downward mass transport, 
respectively, have been set to be /„ = fd{= f ) = 0.5. When this analysis is pursued in the 
framework of a slow mixing process, the physical mechanism of transport and its velocity 
do not play any crucial role; the nucleosynthesis results are then controlled by the path 
integral of reaction rates over the trajectory. If needed, the mixing velocity v can then 
be simply derived as: v = 4 x ■ As discussed elsewhere (Guandalini et al. 2009), 
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very different scenarios may open up when the transport is operated by fast mechanisms, 
suitable to yield mixing time scales comparable to (or shorter than) the ^Be half-life. The 
complex processes associated to the nucleosynthesis of low mass nuclei like D, ^He, Li, Be 
and B will however be postponed to another work. 

After core He-exhaustion, our low-mass stellar models evolve to the AGB stages. 
General discussions of the physics of AGB stars can be found in Iben & Renzini (1983); 
Busso et al. (1999); Straniero et al. (2006). On the AGB, deep mixing phenomena can 
occur again, and there is actually clear observational evidence of their presence in those 
final stages (Busso et al. 2010). In fact, below the convective envelopes of these evolved red 
giants, conditions of chemical homogeneity similar to those found after FDU are estabhshed. 
During the early phase of the AGB (before thermal pulses start) this is guaranteed by 
the previous advancement of the H-burning shell during the first ascent to the RGB. 
Subsequently, when thermal instabilities from the He shell begin, the downward envelope 
expansions (which follow them) easily reach down to the H-He discontinuity, even before the 
penetration into the He-rich layers occurs, a fact that formally defines the third dredge up, 
or TDU. In most of the interpulse periods of AGB phases, therefore, H burning proceeds 
in an initially-homogeneous environment, i.e. in conditions suitable for the occurrence of 
extra-mixing. 

Various experimental constraints (and in particular the spread in the C, O and Al 
isotopic ratios shown by presolar grains) suggest that the details of the transport vary 
from star to star (NoUett et al. 2003). One has also to mention that not all the physical 
mechanisms suitable to drive chemical transport in first-ascent red giants are necessarily 
available also during the AGB stages. This was discussed by Busso et al. (2010), who 
showed that the operation of deep mixing on the RGB, by whatever cause, consumes '^He, 
so that the molecular weight inversion generated by its burning occurs in AGB phases 
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only if the amount of ^He remaining in the envelope is substantial (as is the case in the 
work by Charbonnel & Lagarde 2010). Otherwise, different physical mechanisms must be 
looked for. In Fig.l we better illustrate this point reporting the variation of the molecular 
weight (//) across the radiative region of a 1.5 and and a 2.0 M© star with solar metaUicity. 
While the inversion on the RGB is appreciable in both the considered models, on the 
AGB it is preserved only if a sufficient supply of ^He remains. In particular, for the AGE, 
different lines refer to different abundances of ^He resulting from the previous RGB phase: 
no destruction, a consumption by factors of 3 or 10. It is evident that the stronger the 
previous ^He consumption is, the weaker is the /i inversion. 



5. Results for the RGB phases 

Table 1 illustrates the main changes in the nitrogen and oxygen isotopic ratios induced 
by FDU, for the sample case of a 1.5 Mq star. In order to allow for comparisons with 
previous results and to disentangle the reasons of the changes, the table is presented for the 
solar metaUicity and for the typical cases Z = 0.01 and Z = 0.02, often used in deep mixing 
discussions. As a comparison, the values previously adopted by NoUett et al. (2003), taken 
from Boothrotd & Sackmann (1999) for the same stellar mass are also shown. 

Table 1 



Nitrogen and Oxygen Isotopes at First Dredge Up. 
(1.5 M0 models with updated reaction rates) 



Case 


15N/14N 


170/160 


180/16, 


Z =0.01 


6.17-10-^ 


1.44-10-3 


1.57-10 


Z =0.014 (©) 


6.27-10-^ 


1.21-10-3 


1.59-10 


Z =0.02 


6.65-10-^ 


9.99-10-^ 


1.62-10 



NoUett et al. (2003) LSl-lO'^ S.IO-IO"^ 1.62-10-3 
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Table 2 



Nitrogen and Oxygen Isotopes at FDU 
(Variation with reaction rates, Z — Zq) 

a): updating only the ^^N(p,7)^^0 



Mass (M0) 


15N/14N 


170/160 


180/160 


1.0 


2.07-10-3 


4.29-10-4 


2.42-10-3 


1.2 


1.51-10-3 


5.44-10-4 


2.19-10-3 


1.5 


1.09-10-3 


1.26-10-3 


2.01-10-3 


1.7 


9.16-10-^ 


2.30-10-3 


1.91-10-3 


2.0 


7.79-10-4 


4.87-10-3 


1.83-10-3 


b) ADEIO & ILIIO rates; new [^^ 


'N/^4n]q 


Mass (Mq) 




170/160 


180/160 


1.0 


1.19-10-3 


3.81-10-4 


1.89-10-3 


1.2 


8.64-10-4 


4.82-10-4 


1.73-10-3 


1.5 


6.27-10-4 


1.11-10-3 


1.58-10-3 


1.7 


5.25-10-4 


2.04-10-3 


1.51-10-3 


2.0 


4.46-10-4 


4.32-10-3 


1.44-10-3 



- 18 - 



Table 1 reveals that the estimates for the surface isotopic ratios of nitrogen and oxygen 
at FDU can change by large amounts with respect to previous findings. In particular, this 
is true for ^^N/^^N and for ^^O/^^O. In order to show the selective relevance of the model 
choices in producing these results, Table 2 shows the same nitrogen and oxygen ratios as 
obtained: a) updating, with respect to NoUett et al. (2003), the physical model parameters, 
the value of the solar metallicity and the ^^N(p,7)^^0 rate; and b) adding also all the other 
new reaction rates by ADEIO and ILIIO, as discussed in section 2, and using the new 
(Jovian) value for the solar system nitrogen isotopic ratio. This is done for the masses 1, 
1.2, 1.5, 1.7, 2 Mq, at solar metallicity. 

A comparison between the two tables makes clear that, while for nitrogen the change is 
mainly due to the choice of the new solar value, the large shifts for ^^O/^^O descend from 
both the new abundances and the new rate for proton captures on ^^N. Indeed, both the 
mentioned upgrades induce small but effective changes in the stellar structure, in particular 
in the temperature profile and in the maximum penetration of the envelope during the 
FDU. These changes are obviously larger for stars where CNO cycling is more efficient: 
our 2 Mq model has ^^O/^^O — 5.14x10"^, which is more than twice the value of 2xl0~^, 
previously used in deep mixing models by NoUett et al. (2003) and Boothrotd & Sackmann 
(1999). In this case, the new adopted solar metallicity increases the ^'^O/^^O by about 
65%, while the new ^^N(p,7)^^0 adds a further 35%. In the upper panel of Fig 2 we report 
the temperature profile at central hydrogen exhaustion characterizing four different 2 Mq 
models, run with different inputs: 

1. standard case (we shall call it ST); it has the new ctmn and Z — Zq; 

2. test case with the new an^ and Z—0.02 (this is the solar metaUicity derived by 
adopting the compilation of Anders & Grevesse 1989); 

3. test case with the old ai4^ and Z — Zq\ 
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4. test case with the old 1714^ and Z=0.02. 

The standard case shows the highest temperature profile, this fact leading to a larger ^^O 
production in the region between 0.3 and 0.8 Mq (the inner core is not affected by the 
FDU, while the external layers have temperatures too low to activate H- burning). In the 
lower panel of Fig. 2 we report the abundances of hydrogen, ^^O and ^^O of our ST case. 
In the same panel we also report the maximum penetration of the convective envelope 
during FDU. Note that, although the envelope penetrates less than in the other tests 
(MLtom=0.302, M2„^,„^=0.290, M3„^,„^=0.298, Mf„,,„^=0.286), the standard case shows the 
largest ^^O/^^O value, because of the strong dependence of this ratio on the temperature 
(see also NoUett et al. 2003, their Fig. 8 in particular). 

A look at Table 2 reveals that the other reaction rate changes of section 2, including 
those directly affecting oxygen isotopes, induce smaller (although non-negligible) effects 
than the ^^N(p,7)^^0 does, in altering the isotopic mix of oxygen at FDU. 

Fig. 3 shows this in graphic form: filled squares represent the isotopic composition 
of oxygen at FDU in our models, while open squares show the previous situation, with 
the old NACRE rates and the old choice for the solar metallicity. The new results 
confine the stellar masses that produced presolar oxide grains within a very narrow range 
(1 < M/Mq < 1.5 - 1.7). 

The abundances seen at the surface of a giant star near the tip of the RGB are then 
known to differ from those at FDU for the effects of deep mixing phenomena, which change 
the envelope composition by linking its material to the hotter, radiative regions above 
the H-burning shell, where proton captures re-establish a varying abundance distribution, 
starting from an initially-homogeneous composition. Fig. 3 also shows the effects of deep 
mixing on the RGB (for a specific choice of the parameters defined in Section 4, namely A = 
0.15 and Mq — 0.3). As the figure shows, presolar grains with a moderate ^^O destruction. 
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mainly of group 1 (Nittler et al. 1997), can be reproduced by this rather extreme processing, 
provided the stellar masses are confined within the range mentioned above. This range is 
the same where most MS-S stars are found. We shall come back on this point in the next 
section. Most grain data, however, and especially those with both ^^O/^^O and ^^O/^^O 
ratios below 0.0012, require even more extended processing than presented in the figure and 
clearly indicate that deep mixing cannot occur only on the RGB (see Busso et al. 2010). 

Typical abundance profiles created by nucleosynthesis in the region below the convective 
envelope after the luminosity bump are shown in the four panels of Fig. 4 taken from our 2 
Mq model with initially solar composition. In regions below A ~ 0.2 (corresponding to the 
vertical dashed line) the ^He+^He reaction starts to dominate over the ^He+^He one (that 
induces the inversion of the molecular weight); this is therefore the innermost position from 
which thermohaline diffusion can start. Notice that some remarkable abundance changes, 
including the production of a critical nucleus like ^^Na, and an enhanced destruction of 
^^O, occur immediately below this zone. If only thermohaline mixing were to take care of 
transporting matter to the envelope on the RGB, then we would confirm the recent findings 
by Charbonnel & Lagarde (2010), according to which these nuclei are only moderately 
affected in RGB stages. If instead deeper mixing mechanisms were at play, this result would 
change drastically (as we have shown for the destruction of ^^O in Fig. 3). 

Indication on the mixing depth can obviously be derived by observational data. We 
must however recall that spectroscopic measurements, especially for isotopes, have very 
large uncertainties, so that their constraints should be used with care. Furthermore, the 
metallicity of Population I red giants spans over a large range, while our models cover only 
a factor-of-two interval. As a consequence, among the several existing compilations of CNO 
abundances for red giants of the galactic disk we can consider only those with metallicities 
in the range 0.01 < Z < 0.02, as the interplay of nucleosynthesis and mixing is extremely 
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sensitive to metallicity variations. The best comparisons for our solar-metallicity cases 
{Z = 0.014) should come from M67, whose composition has always been confirmed to be 
very close to solar (Randich et al. 2006). For this cluster CN abundances and ^^C/^^C 
ratios were determined by Brown (1987) and Gilroy & Brown (1991), but the oxygen 
abundance was not measured and rather assumed from older works. New determinations 
of the oxygen concentration in M67 now exist (Tautvaisiene et al. 2000; Yong et al. 2005; 
Pancino et al. 2010) and their [0/Fe] ratios essentially agree with one another; however, 
other rather different suggestions are also available (Priel et al. 2010). We shall adopt here 
the [0/Fe] and [Fe/H] values by Yong et al. (2005), rescaling oxygen for our updated choice 
of the solar abundance from Asplund et al. (2009). The results of our Z = 0.02 case can 
instead be compared to observations in evolved RGB stars of open clusters younger than 
the Sun, hke e.g. those studied by Smiljanic et al. (2009) and Mikolaitis et al. (2010). 
Only few stars from samples of the old disk (see e.g. Cottrell & Sneden 1986; Shetrone et 
al. 1993) can instead be used for our purposes, namely those whose metallicity is close to 
solar. Fig. 5 shows a comparison of our models with observations selected according to 
the above criteria. Given the uncertainty of spectroscopic measurements (especially the 
older ones), the agreement is quite good, especially for M67. The observations do require 
deep mixing at a moderate temperature (A = 0.2) and at circulation rates 0.1 < Mg < 
0.3. From the value of A we infer that thermohaline diffusion might be adequate for this 
purpose. However, warnings on the mixing velocity and on the diffusion coefficient deduced 
from recent literature have to be considered and we refer to section 7 for a more detailed 
discussion of that. 
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6. Results for the AGB phases 

For the late evolutionary stages, the adoption of low temperature C-enhanced opacity 
tables, the choice of a new mass loss law (limiting the number of thermal pulses) and the 
different treatment of the inner convective border all play important roles in controlling 
the final results. While, again, we direct the reader to Cristallo et al. (2009a) for details, 
we need to recall here that the present AGB models have now a much stronger support 
from observations. In particular, the new opacities and the more effective mass loss law 
adopted combine in inducing lower values for the surface temperature and luminosity. In 
their turn, these facts affect the colors, the infrared emission and the luminosity function of 
AGB stars. The new models now succeed in reproducing the photometric and spectroscopic 
properties of their observational counterparts; in particular, the low effective temperatures 
of C stars at different metallicities (Dommguez et al. 2004; Busso et al. 2007b) can now be 
naturally reproduced. The stronger mass loss rates reduce the number of thermal pulses, 
but the introduction of an exponentially decaying profile of convective velocities at the base 
of the envelope increases the TDU efficiency, so that the final enhancement factor of the 
most significant nuclei produced by AGB stars is roughly the same as before; but now this 
is obtained together with a good reproduction of the luminosity functions of AGB stars 
(Cristallo et al. 2011, in preparation). 

We also recall that the adopted algorithm for the velocity profile accounts, after each 
TDU episode, for the formation of a tiny region enriched in protons where, subsequently, 
^^C can be formed. This is the radiative ^^C pocket where {a,n) captures on ^^C provide the 
main neutron source for the synthesis of heavy n-rich nuclei in AGB stars (see e.g. Busso et 
al. 1999; Lugaro et al. 2003; Straniero et al. 2006; Sneden et al. 2008), thus complementing 
the production of lighter n-rich nuclei in massive stars (see e.g. Busso & Gallino 1985; 
Raiteri et al. 1993; Kappeler et al. 1994; Pignatari et al. 2010). We however postpone the 
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analysis of this important point to a future paper. 

Deep mixing in thermally-pulsing AGB stages operates on the composition established 
above the H-burning shell and homogenized after each pulse by the penetration of the 
convective envelope. Typical abundances established in those layers by proton captures 
are shown in Fig. 6. Previous work (NoUett et al. 2003; Busso et al. 2010) leads us to 
expect that extra mixing mainly operates down to regions where A = 0.1 to 0.2. These 
zones can be identified in Fig. 6 from the two vertical lines. As is clear from the plots, 
remarkable nucleosynthesis changes are induced there by the adoption of the new reaction 
rates (continuous curves). Among the differences, we find a larger destruction of ^He, ^^N, 
i8q^ i9p g^jj^ ^ lower production of ^^Na. This last result is directly linked to a lower 
destruction of ^^Ne and ^^Ne. The important radioactive nucleus ^^Al is more effectively 
produced than before, but only marginally. In general, these results can be interpreted 
as due to the slight increase in the temperature induced by the new (less efficient) rate 
for ^^N(p,7)^^0, superimposed to the specific effects of the rate changes for heavier nuclei, 
which have a more local and limited effect. 

As a consequence of the model changes and of the reaction rate updates, the combined 
operation of mixing and nucleosynthesis leads now to abundance changes in the envelope 
like those summarized in Fig. 7 (for the 1.5 Mq model) and in Fig. 8 (for the 2 Mq model). 
For these examples we adopted, on the RGB, A = 0.2, Mq = 0.3; on the AGB, A = 0.1, 
Mq = 1. These choices only serve for the purpose of illustration. A wider sample of cases 
for the AGB is presented in the subsequent Fig. 9 and Fig. 10^. 

The above figures referring to the AGB clearly show a saw-tooth trend in the 

■^For researchers interested to comparisons with specific observations, looking for an even 
wider series of cases, we can provide the nucleosynthesis results upon request. 
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abundances of several nuclei. This is induced by the sudden mixing with He-rich layers at 
the third dredge up (TDU), after each pulse. This trend is therefore pronounced for those 
nuclei that are largely affected by He-shell nucleosynthesis, in particular ^^C, ^^Ne and ^^Al. 

The ^^C-enrichment by TDU is partly compensated by its destruction in deep-mixing 
processing during the inter-pulse stages. The operation of such mixing at even moderate 
efficiencies has the consequence of preventing the formation of a carbon star in our models 
up to 1.7 M0. Combining this with the changes already discussed for the FDU (see Fig. 3), 
wc can say that both MS, S giants and oxide presolar grains should be produced mainly 
by stars below this mass limit. Obviously, a more massive star (2 Mq or more), before 
becoming C-rich, also passes through an 0-rich phase in which TDU induces MS-S spectral 
characteristics; but the high efficiency of dredge-up we have now implies that this is a 
not a dominant phase. This can be inferred from the limited amount of mass lost in the 
0-rich stages by a 2 M0 model star becoming C-rich. For moderate extra-mixing {Mq — 
0.1) this amounts to 0.4 and 0.2 Mq at solar metaUicity and a,t Z — 0.01, respectively. 
By contrast, more than 1 Mq is lost in the C-rich phase. On the contrary, a 1.5 — 1.7 
Mq star looses its whole envelope (at least 0.8 Mq) in the form of 0-rich material. If one 
also considers that the initial mass function favors these lower masses, one can deduce that 
most of the observed MS and S giants should have a mass smaller than for C-stars and 
that among these low mass AGB stars one should find the typical progenitors of presolar 
oxide grains. There is obviously a possibility of finding MS-S giants among intermediate 
mass stars (M > 4M0), where hot bottom burning prevents the formation of a carbon star. 
However, for them the expected isotopic ratios (especially for Oxygen) arc very peculiar 
(e.g. with large excesses of ^^O). Iliadis et al. (2008), by comparisons of intermediate-mass 
models with presolar grain abundances, convincingly excluded any connection of IM stars 
with presolar grains; presently, we don't know of even one single grain really requiring an 
intermediate-mass progenitor. 
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We can see this in some more detail. As already stated in discussing the RGB phases 
and in Busso et al. (2010), oxygen isotopes in presolar oxide grains provide the most 
stringent constraints to stellar models with deep mixing, due to the fragility of ^^0 and 
to the sensitivity of ^^O to temperature. Fig. 11 (left panel) presents these constraints, 
showing the oxygen isotopic ratios in group 1 and group 2 grains together with model 
curves. These last include typical trends for extra-mixing on the RGB at a moderate 
efficiency, plus the results of deep-mixing computations for the AGB, for various choices 
of the parameters. Notice that, by including the 1 M© star, the whole range of ^'^O/^^O 
ratios is covered. This is another effect of the model changes and corrects previous findings, 
based on older models, according to which the grains at the extreme left would not be 
explained by deep mixing at solar metallicity (Wasserburg et al. 1994; NoUett et al. 2003). 
The main reason of the change lays in the more marked destruction of ^^O at moderate 
temperatures, where ^^O is still untouched (see Fig. 4, second panel). This is found also 
in very low-mass models and allows their deep-mixing curves to proceed almost vertically, 
covering a previously unaccounted area (not only the new rate for the ^^N(p,7)^^0, but also 
those for proton captures on ^^O are important in determining this result). The left panel 
of Fig. 11 also shows that the whole set of abundances for ^^0-poor grains (with ^®0/^^0 
ratios below 0.0015) is explained by models for stars below about 1.7 Mq, those that will 
never become G rich. 

The right panel of Fig. 11 then presents the ^^O/^^O ratio as a function of ^^Al/^^Al, 
for those grains that show ^^Mg excesses from ^^Al decay. Here Mq is assumed to be 3 (a 
rather efficient case). Again, most of the experimental results can be accounted for. A 
small group of grains, for ^''O/^^O ratios below 0.001 and ^^Al/^^Al above 0.01, are however 
not explained. We notice here that no modification to the model might cover this area. In 
fact (as an example) mixing at higher temperature would improve the fit to ^^Al but would 
worsen the problem for ^^O (not to mention the fact that mixing at higher temperatures 
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would affect the stellar structure). There appears to be a unique way out: an increase 
in the reaction rate for the ^^Mg(p,7)^^Al reaction at temperatures sufficiently small to 
minimize the effects on ^^O. In the energy region of our interest the main contribution to 
that reaction rate is due to resonances (at 58 and 93keV) for which large uncertainties are 
still present (Ihadis et al. 1996). We found that, by artificially increasing their strengths 
by a factor of 5 with respect to present recommendations, grains rich in ^^Al and not in 
^''O can be accounted for (e.g. using A — 0.05 and Mq = 3 in a 1.2 Mq stellar model). We 
leave this suggestion for future verifications. Recent new measurements (especially from 
LUNA) are still under analysis, so that we should not wait a long time before knowing if 
this prediction can be confirmed. 

It is important to notice that the same presence of ^^Al in the grains, at a level higher 
than allowed by TDU (which therefore requires the operation of deep mixing), poses crucial 
constraints to the physical mechanisms driving the transport of processed matter. Contrary 
to what was found on the RGB, here thermohaline diffusion cannot be an explanation. Not 
only its effectiveness is put in doubt (or at least reduced) by the previous ^He consumption, 
but the layers where it starts are too cool to host any ^^Al. This is evident in Fig. 6. On 
the AGB, the ^He+^He reaction prevails over the ^He+^He for A > 0.15, identified in the 
first panel of the figure by the kink in the ^He curve. As the third panel shows, the whole 
abundance of ^^Al is produced inside that zone. As we showed in Fig. 11, accounting 
for ^^Al in oxide grains requires A < 0.1. We believe this is a clear demonstration that 
thermohaline mixing is not the mechanism driving extra-mixing in AGB stars. Even its 
action as a driver of magnetic-like mechanisms (Denissenkov et al. 2009) would not be useful 
in this case. Completely different physical processes, like e.g. magnetic buoyancy (Busso et 
al. 2010) or gravitational waves (Denissenkov & Tout 2003), must necessarily be at play. 

Further comparisons with observations are provided by Fig. 12. Its left panel shows 
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C/0 and ^^C/^^C measurements in both 0-rich (Smith & Lambert 1990) and C-rich stars 
(Abia et al. 2002, 2010). The curves refer to our 1.5 and 2 Mq models; the more massive 
case is the one that reaches C/0 > 1; its C-rich phases are marked by a continuous hne, 
the 0-rich ones by a dashed hne. We adopted A = 0.2 and 0.25 and Mq values as indicated 
by the labels. It is clear that the whole area of the observations is very well accounted for. 
This figure offers a rather global test from stars, complementing those derived from oxide 
grains; it confirms that the general scenario is correct and that MS, S and C stars have 
indeed the masses we suggested and do experience extra-mixing. 

A much less satisfactory comparison is presented in the right panel of Fig. 12, where 
our model curves are superimposed to the the nitrogen and carbon isotopic ratios measured 
in SiC grains of AGB origin (see e.g. Zinner 2005; Zinner et al. 2006, and references therein). 
The models in the figure are labeled by the values of A and Mq. They mainly cover the area 
of the so-called "mainstream" SiC grains (with carbon isotopic ratios between 20 and 100), 
but only for the upper part of the range in ^^N/-'^^N. Models with rather extreme Mq values 
can reach the area where the grains with the lowest ^^C/^^C ratios (called A & B grains) 
are present, but again only for very high values of ^^N/^^N. These inconsistencies are not 
new and have characterized the attempts at explaining the nitrogen isotopes in SiC grains 
for many years (see e.g. NoUett et al. 2003). Notice that the two filled circles, at ^^N/-^^N = 
680 and 5300, refer to the only (very uncertain) data from red giants (Wannier et al. 1991) 
and confirm that the RGB ratios must indeed be much larger than solar. The low values of 
many SiC grains have no explanation at all in low mass stars. In principle, in the so-called 
^^C-pocket of AGB stars a reservoir of ^^N exists. This is generated trough proton and 
neutron captures on ^^N and later, at the onset of the thermal pulse, by the activation of 
the ^^0(p,q;) reaction, see also Abia et al. (2010). However, these zones are not expected to 
be ejected and their ^^N is subsequently destroyed in the convective thermal pulse. 
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Contamination of the grains with solar-system nitrogen seems unhkely, as there is no 
correlation between the ^^N concentration and grain size, as would be expected (Zinner, 
private communication). A possible alternative is the contamination from cosmic ray 
spallation, where ^^N is produced abundantly. In such a case the ^^N content would be a 
measure of the grain age and would not inform on stellar nucleosynthesis. 

7. TowEird a physical model of deep mixing: constraints on the parameters 

As a concluding remark we want to make some comments on the recent numerical 
simulations of extra-mixing presented by Denissenkov (2010) and by Denissenkov & 
Merryfield (2010). The work was performed through an analogy with current research on 
thermohaline processes known to occur in the oceans. The models were then extrapolated 
to the stellar case, on the basis of the values of physical parameters valid for red giants; 
this treatment has obviously a much higher degree of sophistication than in simple ID 
analysis, but still includes very uncertain parameters. Anyhow, the authors suggest 
that the molecular weight inversion generated by ^He burning induces finger-like mixing 
structures whose average aspect ratio (length over diameter, a = I /a) is too small by a 
factor of at least 50 as compared to what observations require; this aspect ratio, squared, is 
proportional to the diffusion coefficient in a diffusive treatment, or to the mixing velocity in 
a circulation- like transport (Nollett ct al. 2003). In Denissenkov (2010) it is also mentioned 
that, for explaining the small number of red giant stars that are very Li-rich (loge(Li) 
> 1.5) the aspect ratios found with thermohahne processes would be too small by a factor of 
several hundreds; in fact, producing ^Li is known to require a kind of transport fast enough 
to save to the envelope the short-lived ^Be (ti/2 = 53 days). These results would indicate 
that thermohaline diffusion alone is insufficient to drive the chemical mixing required in red 
giants; as a way out, in Denissenkov & Merryfield (2010) it is suggested that a couphng 
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with magnetically-induced effects might be more promising. 

The above results should now be verified by other authors, before final conclusions 
can be drawn. Recent works either attribute the isotopic mix of red giants to the effects 
of thermohaline mixing alone (Charbonnel & Zahn 2007) or confirm that models for this 
process find velocities that are too low by a large factor, thus opening the possibility that 
other mechanisms contribute (Cantiello & Langer 2010). 

Although our simple parameterized approach cannot enter into such details, we found 
the treatment by Denissenkov (2010) and Denissenkov & Merryfield (2010) very instructive 
as a tool to infer the required properties of suitable physical models for extra-mixing. In 
particular, it was suggested (Busso et al. 2007a; Denissenkov et al. 2009) that magnetic 
buoyancy might offer an alternative to thermohaline diffusion, with the advantage that it 
would not require a nuclear process as a trigger, hence it would have a wider possibility 
of application (e.g. to the formation of the ^^C-pocket). If magnetic buoyancy were to 
proceed through the development of 0-shaped loops, formed as instabilities in MHD waves 
(Parker 1994), then each such structure would look hke a couple of mixing fingers. Let us 
make two illustrative examples to derive limits on the aspect ratios of these structures. As 
a zero-order approximation we can consider the model by Busso et al. (2007a), neglecting 
any slowing process (such as phase mixing or thermal exchanges, see e.g. Spruit 1999, 
2002; Denissenkov et al. 2009). In this case the buoyancy velocity is close to the Alfven 
velocity and the aspect ratio, for the cases discussed by Busso et al. (2007a), turns out 
to be in the range 200—300, suitable for Li production. Alternatively, a less efficient 
transport is obtained if the instabilities can travel only a limited distance I before interfering 
destructively (by phase mixing) with previously-detached buoyant structures. The length 
covered is / = TrAt^/t, where t is the available time for the mixing and At a is the crossing 
time of Alfven waves over the star's dimension (Spruit 1999; Busso et al. 2004). In such a 
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case, using the parameters of RGB phases by Busso et al. (2007a), one can easily derive 
that the aspect ratio becomes very small (of the order of 0.1—0.2), suitable to account for 
the requirements of the carbon isotopic ratios in red giants and with a correspondingly 
small velocity. A similarly small velocity is obtained for the magneto-thermohaline mixing 
by Denissenkov et al. (2009). These are only very simphfied, order-of- magnitude estimates, 
derived just by applying to magnetic buoyancy some of the considerations by Denissenkov 
(2010). However, they are already sufficient to underline that this kind of mechanism 
deserves now to be carefully modeled, as it seems to be able to cover the whole range 
of mixing velocities (or aspect ratios) required by virtually all the chemical peculiarities 
displayed by red giant stars. 



8. Conclusions 

In this paper we presented revised parametric calculations of nucleosynthesis in 
extra-mixing phenomena occurring in low mass evolved stars, from the phase of the FDU 
up to the end of the AGB stage. The calculations were based on a recent revision of the 
FRANEC code and on a new extra-mixing code. They incorporated several physical and 
numerical improvements, among which new opacities for carbon-rich envelopes, a new 
algorithm for the determination of the convective border position, and a new mass-loss law. 

With these tools we explored the coupling of nucleosynthesis and non-convective mixing 
in RGB and AGB phases, comparing the results obtained with the old NACRE compilation 
of reaction rates with those descending from the adoption of new recommendations from 
ADEIO and ILIIO. In particular, the specific choices made for each updated reaction rate 
were commented extensively in the Appendix. 

It was found that, among the new reaction rates, the ones inducing the largest 
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variations in the evolution of LM stars are those for proton captures on ^^N and (to a 
lesser extent) for proton captures on -'^'''O. Similar conclusions were recently obtained in 
another context by Magic et al. (2010). We obtained huge changes in the oxygen isotopic 
composition of the envelope, starting from FDU. Their effects include a drastic shrinking 
of the mass interval in which AGB stars can produce oxide grains like those recovered in 
pristine meteorites. Most such grains appear to come from stars below 1.5 Mq. Another 
important revision of previous knowledge was the discovery that all ^*0-poor grains 
(those with ^^O/^^O below 0.0015) can now be explained, without the "forbidden regions" 
previously found for the lowest ^'^O/^^O ratios. Concerning the evolution of AGB stars of 
low mass and solar metallicity we suggested that most MS, S stars should have masses 
below 1.5 — 1.7 Mq] in this mass interval they would never produce C-stars, also thanks to 
^^C processing in deep mixing. 

Constraints coming from observations of CNO nuclei in both RGB and AGB stars were 
well reproduced by deep mixing models with the new rates. In particular, for RGB stars we 
found that extra-mixing processing does not need to penetrate to very high temperatures, so 
that thermohaline diffusion might be an adequate scenario for those phases, if the problems 
recently mentioned in the literature for its too low diffusion coefficients can be solved. We 
then confirmed that deep mixing is required also on the AGB. In this case, the abundances 
of ^^Al in presolar grains indicate that the transport must reach much hotter layers than in 
RGB stages. This fact (and the previous consumption of ^He) induced us to conclude that 
thermohaline diffusion is not a viable explanation for the transport processes occurring in 
the latest evolutionary stages. The existence of a small population of oxide grains with low 
i7q^i6q g^j^j high ^^Al/^^Al ratios was suggested to require some further revision of the 
rates for ^^Al synthesis at low temperatures, probably an increase for the ^^Mg(p,7)^^Al 
rate. 
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The longstanding impossibility of accounting for the whole range of the nitrogen 
isotopic ratios in SiC grains remains unchanged after our work. It was argued that no 
stellar process occurring in low mass stars can solve this problem and that SiC grains might 
be contaminated by cosmic ray spallation processes in the interstellar medium. 

We are indebted to many colleagues for helping us in the understanding of the subtleties 
of low mass star nucleosynthesis. While recognizing the role of all of them is impossible, 
we cannot avoid acknowledging the fundamental lessons we learned from friends like G.J. 
Wasserburg, R. Gallino, O. Straniero and K.M. NoUett. M.L. acknowledges the support of 
the Italian Ministry of Education, University and Research under Grant No. RBFR082838 
(FIRB2008) 

A. The choice of nuclear reaction rates 

A.l. The reactions ^^N(p, 7)1^0 and ^^N{p,^y^O. 

The ^^N(p, 7)^^0 reaction is the slowest reaction of the CN cycle, hence the one 
determining the duration of the process. The low-energy cross section is dominated by a 
resonance at 259 keV. In particular, at the temperatures of interest here (Tg < 0.06), the 
reaction rate mainly depends on the resonance tail, the Gamow window extending (at most) 
to around 110 keV. Such a region is only partly covered by experimental data, spanning 
energies as low as 70 kcV, so the S-factor is very sensitive to the extrapolating curve and 
can be severely affected by the —504 keV subthreshold resonance. In the last years, direct 
measurements (in both surface and underground laboratories) and indirect measurements 
(Dopplcr shift attenuation. Coulomb excitation. Asymptotic Normalization Constant) 
have been performed, primarily triggered by the discovery that the contribution of the 
ground-state transition had been previously overestimated. Prom these measurements. 
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ADEIO recommend S{0) = 1.66 ± 0.12 keVb, which is about half the NACRE value 
-S'(O) = 3.2 ± 0.8 keVb. This entails a significantly lower reaction rate in the range of our 
interest, as shown in Fig. 13, left panel. 

The ^^N(p, 7)^^0 reaction fixes the leak rate from the CN-cycle to the CNO bi-cycle. 
Besides the direct capture contribution, the cross section is characterized by the interference 
of two 1~ resonances at 312 and 964 keV. As for ^^N(p, 7)^^0, the largest energy of interest 
is ~ 110 keV, thus the extrapolation procedure is critical. In NACRE, a very large S(0) was 
obtained (64 ± 6 keVb) while ADEIO recommend 36 ± 6 kcVb following the results from 
an underground measurement (which does not cover the Gamow window anyway) and an 
ANC measurement, yielding an absolute value of the direct-capture contribution ten times 
smaller than NACRE. Consequently, a 30-40% smaller reaction rate is now obtained in the 
temperature region we are dealing with (Fig. 13, right panel), by using the low-energy 
S(E)-factor recommended by ADEIO. 

A.2. The reactions i60(p,7)i^F, i^0(p,7)i«F, and i^O(p,Q;)i%. 

The operation of the ^^0(p, 7)^^F reaction contributes to estabhsh the ratios of oxygen 
isotopes in red giants and plays therefore a major role in the analysis of oxygen-rich presolar 
grains of circumstellar origin. At the energies of interest for this work, and up to ~ 120 
keV, the cross section is determined by the direct capture to the ground and first excited 
states of ^^F. Direct measurements exist down to about 100 keV, so extrapolations are 
needed. In the previous evaluation by NACRE, the greatest part of the error was attributed 
to the model assumed in the calculation. Improved calculations have now reduced the 
uncertainty affecting the S(0) parameter, from 30% to 7.5% (see ADEIO for details), while 
the recommended value is in substantial agreement with the one from NACRE. Fig. 14 
(left panel) shows the resulting reaction rate. 
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Much like the previously- mentioned reaction, also ^^0(p, 7)^^F is expected to modify 
the relative abundances of oxygen isotopes. In the Ecm ^120 keV energy range, the cross 
section shows a resonance at 65.1 keV that strongly influences the reaction rate above 
Tg — 0.02. At lower temperatures, the dominant role is played by direct captures, while 
the subthreshold resonance at -3.1 keV does not provide any significant contribution. The 
65.1 keV resonance and the direct capture term have been the subject of several improved 
calculations, leading to a change in the resonance energy (it was 66 keV in NACRE) and 
in its strength (which is a factor ~ 3.7 smaller than given in NACRE). Concerning the 
direct-capture contributions, the revised calculations yield a ~ 30% smaller S(E) factor 
below 120 keV, which is responsible for the reduced reaction rate below Tg = 0.02, as shown 
in Fig. 14, central panel. 

The ^^0(p, q;)^^N competes with the ^^0(p, 7)^^F channel, removing ^^O nuclei and 
turning the nucleosynthesis path back towards nitrogen. In contrast with the radiative 
capture channel, the —3.1 keV subthreshold resonance gives a major contribution to the 
reaction rate for Tg < 0.02. Updated calculations for these resonances are given in both the 
reviews we are considering. In particular, the -3.1 keV-state proton width is lowered by a 
factor ~ 3 as compared to NACRE. No expression for the reaction rate is given in ADEIO; 
hence we use the reaction rate table in ILIIO, extrapolating down to the Tg = 0.004 — 0.01 
temperature range. A remarkably lower reaction rate is thus obtained (see Fig. 14, 
right panel). Regarding the ^^0(p,Q;)^^N reaction, despite the extensive experimental 
investigation in the last years (ADE10,ILI10, La Cognata et al. 2008, 2010b) has produced 
a significant change in the recommended reaction rate (about 30%) and a reduction of its 
uncertainty (about a factor of 8) at low temperatures, its effect on AGB nucleosynthesis 
turns out to be negligible. The updated reaction rate has been included anyway in the 
present work. 
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A.3. The 2°Ne(p, -ff^Na, ^^Ne{p, -ff^Nei and ^^Ne{p, -ff^Nei reactions 

The ^°Ne(p, 7)^^Na reaction is the starting point of the NeNa cycle and plays a key 
role, being the slowest reaction of the cycle. Inside the Gamow window, extending up to 
~ 130 keV, the astrophysical factor is dominated by a subthreshold resonance laying at 
about -7 keV. Additional 8 resonances show up at higher energies (up to about 2 MeV). 
In ADEIO both the resonance excitation energy and its strength have been re-evaluated 
using more recent measurements, resulting in a small increase of the S(E)-factor, by 2.9%. 
Moreover, a new fit of the energy trend for the S(E)-factor has been provided for the 
numerical integration of the reaction rate. The resulting rate is displayed in Fig. 15, left 
panel, which shows a substantial agreement with the NACRE one and a shghtly reduced 
uncertainty. 

For the '^'^Ne{p,^y'^Na reaction, we adopt the recommendation by ILIIO. They 
make use of nuclear-physics input parameters from Iliadis et al. (2001). Resonances were 
measured for Ecm > 121 keV, corresponding to the upper edge of the Gamow window. Two 
threshold states, corresponding to resonance energies of Ecm — 17 and 95 keV, have been 
considered. While the latter (for which only an upper limit of the spectroscopic factor is 
known) adds a negligible uncertainty to the reaction rate, the former can influence it at 
low temperatures (below Tg ~ 0.02). Though the reaction-rate accuracy has been improved 
with respect to NACRE around Tg ~ 0.03, for Tg < 0.02 the total rate remains uncertain 
by several orders of magnitude (see Fig. 15, central panel). 

The rate of the ^^Ne(p, 7)^^Na reaction is taken from ADEIO and is shown in Fig. 15, 
right panel. The main difference with respect to previous compilations is the treatment 
of the threshold states. Concerning the 151-keV resonance, its spectroscopic factor has 
been considered as a mean value rather than an upper limit. Lower energy resonances, 
whose existence had been considered as tentative, were disregarded by ADEIO. As a result. 
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the recommended rate is smaller than in NACRE above Tg ~ 0.03, but remains in good 
agreement with the previous indications elsewhere. The error has been instead significantly 
reduced. 

A.4. The 23Na(p,7)^^Mg and '^^Na{p,af°'Ne reactions 

The Gamow window of the ^^Na(p, 7)^^Mg reaction extends up to about 140 keV. 
Therefore, the main contribution to the low-temperature reaction rate is due to the 
threshold states. ILIIO mostly make use of the updated nuclear physics input parameters 
by Hale et al. (2004) in this energy region, which have been indirectly measured by means 
of (^He, d) transfer reactions. The contribution of direct capture to the cross section has 
been taken from the same work. The resulting reaction rate is compared with the NACRE 
choice in Fig. 16 (left panel). The updated value lays well within the upper and lower limits 
set by NACRE, although the uncertainty is greatly reduced. 

Concerning the ^^Na(p, a)^'^Mg reaction, the same considerations as for the 
^^Na(p, 7)^°Ne can be repeated. The resulting reaction rate is given in Fig. 16, right panel. 
The new rate is about 100 times smaller than the one recommended by NACRE at the 
lowest temperatures, while the uncertainty is significantly smaller. 

A.5. The reactions ^^Mg{p,-ff^A\, ^'^Mg{p,-f)'^'^Algs, '^'^Mg{p,-ff^A\^, ^^Mgip^-ff^ Al. 

The ^"^Mgd?, 7)^^A1 reaction influences the nucleosynthesis of magnesium and aluminum, 
which are key elements in the analysis of ancient meteorites and presolar grains. No 
experimental data are available inside the Gamow window (which extends up to Ecm ~ 150 
keV). The most recent measurements provide the cross section down to 200 keV (see 
ILIIO and reference therein), where 9 resonances occur (up to 2.3 MeV). In the review by 
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ILIIO, the reaction rate is calculated from the these data, including re-normalizations of 
the strengths of higher-energy resonances and updated Q-values. Moreover, with respect 
to the NACRE compilation, a more reliable direct capture contribution to the S(E) factor 
has been estimated, as well as the contribution of the 214 and 402 keV resonance tails, by 
means of numerical integrations. Despite all these improvements, the resulting reaction 
rate (Fig. 17, left panel) is in overall fair agreement with the NACRE recommendations. 
It is shown for Tg > 0.01, but the reaction rate is there ~ 10~^^cm^mol~^s~^, so that the 
^^Mg(p, 7)^^A1 reaction does not contribute to nucleosynthesis at very low energies. 

As for the previous case, the ^^Mg(p, 7)^^Alps reaction (where "gs" stands for "Ground 
State" ) is important for us to understand the magnesium and aluminum nucleosynthesis in 
evolved red giants. Data exist in a wide energy range, starting from 37 keV and extensively 
overlapping with the Gamow peak. Concerning the nuclear physics inputs, the differences 
of the ILIIO review with respect to NACRE are: (i) the use of re-normahzed strengths for 
resonances above Ecm — 189 keV; (ii) the adoption of the more accurate ground-state 7-ray 
branching ratios /o for the 189, 244 and 292 keV resonances; and (iii) a new calculation of 
the contribution of the 37 keV resonance. The reaction rate is displayed in Fig. 17, right 
panel: the recommended value is in good agreement with NACRE, but the uncertainty has 
been reduced by more than a factor of 2 at the lowest temperatures that are most important 
for us. 

The same considerations done above apply to the ^^Mg(p, 7)^^Alf„ reaction (where "m" 
indicates the isomeric meta-stable state), but in addition the isomeric-state 7-ray branching 
ratio 1 — /o has been taken into account. Once again, while the recommended value by 
ILIIO is in overall agreement with NACRE, the uncertainties are reduced by more than a 
factor of 2 at the temperatures of AGB stars (Fig. 18, left panel) 

For the ^^Mg(p, 7)^'''A1 reaction. Fig. 18 (right panel) shows a remarkable difference 
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between the updated reaction rate by ILIIO and the NACRE value, ranging from a factor 
of 2 up to 8, while the error has been reduced, especially at high temperatures. Indeed, a 
different direct-capture component was adopted, the strengths of the measured resonances 
(above Ecm ~ 100 keV) were re- normalized and the contributions of the lower-energy 
resonances at Ecm = 53 keV and Ecm = 105 keV were re- analyzed. One has to remark that 
some ambiguity on the assignment of quantum numbers for these levels still exists (see 
ILIIO for a more detailed discussion). 

A.6. The reactions 26Al(p, 7)27si, ^'^ A\{p, 7)2831, and ^'^ k\{p, af ^Mg. 

The 2^Al(p, 7)2^Si reaction is important for establishing the equilibrium abundance of 
the radioactive isotope ^^Al, playing a major role in presolar grain studies, 7-ray astronomy 
and solar system formation (see e.g. Busso 2010; Busso et al. 2003; Wasserburg et al. 2006, 
for this last point). Fig. 19 (left panel) shows that the updated reaction rate by ILIIO lays 
within the fiducial interval by NACRE, but at low temperature the recommended value is 
strongly decreased. Also the uncertainty has been greatly reduced. In the estimate of the 
reaction rate, ILIIO introduced a calculated direct-capture component and the contribution 
of 19 resonances. In particular, for almost all resonances below Ecm ~ 500 keV, the 
resonance energies were computed from updated excitation-energy determinations, while 
improved calculations of the proton widths of unobserved low-energy resonances were used 
to evaluate their contribution to the rate. 

For the 2''Al(p, 7)^^81 reaction, ILIIO made use of the nuclear physics input parameters 
reviewed in Iliadis et al. (2001). With respect to the NACRE compilation, the direct-capture 
contribution was recalculated, the strengths of the measured resonances were re-normalized 
(above Ecm ~ 200 keV, to be compared with the upper limit of the Gamow window, 
~ 160 keV) and the contribution of four sub-threshold resonances were re-evaluated (the 
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most important are at Ecm — 72 keV and Ecm — 85 keV). In particular, an improved 
determination of their spectroscopic factors was performed, based on the re-analysis of the 
^^Al(^He, (i)^*Si stripping data. Consequently, the recommended total rate deviates from 
NACRE up to a factor of 20 in the energy range common to AGB stars (Fig. 19, central 
panel), mainly because of improved estimates of spectroscopic factors for the threshold 
states. 

A total of 91 resonances at energies of Ecm — 72 — 2966 keV were taken into account 
by ILIIO for calculating the reaction rate for the process ^''A^p, a)'^^Mg. For the measured 
resonances (for E^m ^ 486 keV), the same nuclear physics input parameters considered in 
Ihadis et al. (2001) were adopted. The same holds for threshold resonances, whose strengths 
were calculated from a re-analysis of the ^'''Al(^He, (i)^^Si stripping cross section and from 
spectroscopic data in the literature (see ILIOl and references therein). The improved 
estimates of spectroscopic factors for the threshold states, as in the ^^Al(p, 7)^^Si case, 
determine an increase of the recommended total rate with respect to NACRE by a factor of 
10 (Fig. 19, right panel). 

Finally, we recall that (p,7) and {p,a) captures on fluorine have not been addressed by 
either ADEIO or ILllO. However, recent measurements (Couture ct al. 2007; Spyrou et al. 
2000) only confirmed previous suggestions, so that the rates do not require updates, and we 
adopt them from NACRE. 
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Fig. 1. — Relative variation of the molecular weight (Afx/fx) for a 1.5 and and a 2.0 Mq star 
with solar metallicity, in the layers hosting ^Hc burning above the H shell. The n inversion, 
present on the RGB, is preserved also on the AGB only if a sufficient supply of ^He remains. 
For the AGB, different hues refer to different abundances of ^He resulting from the previous 
RGB phase, as indicated by the labels. 
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Fig. 2. — Upper panel: the profile of the stellar temperature as a function of the fractional 
mass at central H exhaustion in a model star of 2 Mq, run with the indicated choices of the 
reaction rate for proton captures on ^'^N and of the solar metallicity. Znew and Zoid refer to 
the values reported by Asplund et al. (2009) and Anders & Grevesse (1989), respectively. 
Lower panel: the ^^O and ^^O abundance profiles at the same moment to which panel 1 
refers. The plot is for our standard case (i.e. has the new ^^N(p,7)^^0 reaction rate and the 
new solar metallicity). The arrow indicates the layer reached by FDU. 
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Fig. 3. — The oxygen isotopic mix in the envelope of our solar-metaUicity stars at FDU (sohd 
hne). Full squares indicate the mass of the model star. As a comparison, open symbols (see 
arrows for associations) show the composition for the masses 1.2, 1.5 and 2 M0 in the old 
scenario with NACRE rates and with the solar metallicity from Anders & Grevesse (1989). 
Dashed hues (and the corresponding shaded area) show the range of values covered by deep 
mixing on the RGB, for A = 0.15 and Mq = 0.3. These are rather extreme values of the 
parameters for the RGB phase, aimed at showing the maximum effects. The grey data 
points refer to measurements in presolar grains (Choi et al. 1998; Nittler et al. 1997, 2008, 
WUSTL Presolar Database http : / /presolar. wustl.edu/ pgd/). We plot those of group 1 
(open circles) and of group 2 (filled circles). The solar symbol shows the initial isotopic mix. 
It is clear that extra-mixing on the RGB, even operating with extreme efficiency, does not 
suffice to explain the data of group 2 grains. 
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Fig. 4. — Distribution of abundances for some light and intermediate-mass isotopes in the 
radiative layers between the H-burning shell and the convective envelope. The plots are 
taken shortly after the moment when shell-H burning reaches the regions homogenized by 
FDU and re-establishes chemical gradients there. Dotted lines in the figure represent the 
old results, obtained before our reaction rate updates, while continuous lines represent our 
new estimates. The vertical dotted line and the dashed one indicate, respectively, the place 
where A = 0.2 and where the energy from the H-burning shell is maximum (Th). (A color 
version of this figure is available in the online journal.) 
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Fig. 5. — Deep-mixing results for CNO nuclei, as compared to observations in RGB stars of 
metallicities close to solar. Panel a): model curves correspond to A = 0.2; they proceed from 
right to left, starting from abundances at FDU. Continuous lines refer to Mq = 0.3, dotted 
lines to Mg = 0.1, dashed lines to Mg = 0.03. Black stars show post-FDU stars from M67 
(Brown 1987; Gilroy & Brown 1991); asterisks are data from Cottrell & Sneden (1986); open 
squares refer to open clusters by Mikolaitis et al. (2010). Circles refer to data from Smiljanic 
et al. (2009), in particular to NGC5822 (filled symbols) and NGC2360 (open symbols). Only 
solar metallicity models are shown, to avoid excessive crowding; the cases at Z — 0.02 would 
be slightly shifted upward. The typical uncertainty (la) is indicated. Panel b): same type 
of comparison, this time for the carbon isotopic ratio versus C/0. Models of Z = 0.02 and 
Z — 0.01 have been included. Panel c): same comparison, for the N/0 versus C/0 ratios. 



- 54 - 



LogT 

7.8 7.7 7.6 7.5 7.4 7.3 



LogT 

7.8 7.7 7.6 7.5 7.4 7.3 



:--io- 



§'-12 
—I 

-14 
-16 
-18 
-20 



1 1 1 1 1 

1 


1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 


■1 




c 






'Be 


-\ 
-\ 
1 








^o,.i 


A=P,.2,,,,,,, 




-5.2 -5.1 -5 -4.9 -4.8 -4.7 -5-2 -5.1 -5 -4.9 -4.8 -4.7 



Log(R/R^ 
LogT 



Log(R/R^ 
LogT 



7.8 7.7 7.6 7.5 7.4 7.3 7.8 7.7 7.6 7.5 7.4 7.3 

-4r 



-3 
-4 
^ -5 

o) -6 

o 

—I 

-7 

-8 
-9 





1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 






































: , IT,H, 1 , 







O -7 




I ' ' ' ' I 



25, 



Mg 



21 



Ne 



A=P.-1A,?Q.,2 



-5.2 -5.1 -5 -4.9 -4.8 -4.7 -5.2 -5.1 -5 -4.9 -4.8 -4.7 
Log(R/RJ Log(R/RJ 



Fig. 6. — Distribution of abundances for some light and intermediate-mass isotopes in the 
radiative layers between the H-burning shell and the convective envelope. The plots are 
taken roughly at mid-AGB evolution. The two vertical dotted lines and the dashed one 
indicate the places where A = 0.2, 0.1 and 0.0, respectively. This last position corresponds 
to the layer where the energy from the H-burning shell is maximum (and T = Th). (A color 
version of this figure is available in the online journal.) 
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Fig. 7. — An example of the evolution in time of the abundances by mole in the envelope 
for representative isotopes. The RGB and AGB phases of a 1.5 star of solar metallicity 
are shown for a specific choice of the mixing parameters. (A color version of this figure is 
available in the online journal.) 
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Fig. 8. — Same as in the previous figure, but for a 2 M© model star. (A color version of this 
figure is available in the online journal.) 
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Fig. 9. — Evolution in time of the abundances by mole in the envelope of a 1.5 Mq star, 
during the AGB stages. The plot considers three values of the mixing parameter Mg and 
two choices of the temperature parameter A. (A color version of this figure is available in 
the online journal.) 
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Fig. 10. — Same as in the previous figure, for a stellar model of 2 M0. (A color version of 
this figure is available in the online journal.) 
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Fig. 11. — The curves show the oxygen and aluminum isotopic mix in the envelopes of our 
solar-metallicity stars during advanced evolutionary stages, as compared to presolar oxide 
grain compositions (WUSTL Presolar Database http : / /presolar. wustl.edu/ pgd/). Left 
panel: FDU is indicated by the heavy, long-dashed curve; short-dashed lines and the shaded 
area show the effects of deep mixing on the RGB, for A = 0.2 and Mg = 0.1. This is a more 
moderate choice than in Fig. 3. The symbols for data points are the same as in Fig. 3. 
Continuous lines refer to model results for extra-mixing on the AGB (all models for a specific 
mass are named a), b), c), d), e) from left to right). They correspond to specific choices of 
the parameters, namely: a) A = 0.1, Me = 1; b) A = 0.1, Mq = 3; c) A = 0.15, Mq = 1; 
d) A = 0.15, Me = 3; e) A = 0.2, Mq = 1. Models for M < I.TMq explain essentially all 
the data for ^^0-poor grains. In particular, group 2 grains mainly derive from very low-mass 
stars (below 1.5 Mq). Right panel: the ^^O/^^O ratio as a function of the ^^Al content. 
Continuous lines refer to A = 0.05; long-dashed lines to A = 0.1; short-dashed lines to A = 
0.15; dotted lines to A = 0.20. Model curves proceed from left to right, converging toward 
the equilibrium ^''O/^^O at the specific temperature of each model. Most of the data points 
are accounted for (see text for the parameter values). 
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Fig. 12. — Left panel: observations of C/0 and ^^C/^^C ratios in AGB stars of types MS, 
S (white triangles), SC (open circles) and C(N) (filled circles), as compared to a few of our 
curves from models of 1.5 (grey curves) and 2 Mq (black curves) at solar metallicity, with 
extra-mixing and TDU. Continuous lines refer to the C-rich phases, dashed hues to the O- 
rich ones. The labels indicate the choices for A and Mq. The area of the data is well covered 
by the models, indicating that the parameter choices should be rather typical of real AGB 
stars. Right panel: The ^^N/^^N ratios of SiC grains recovered from pristine meteorites, as a 
function of their ^^C/^^C ratios. Open symbols represent the A+B grains, grey symbols the 
so-called "mainstream" ones (Zinner 2005; Zinner et al. 2006, WUSTL Presolar Database 
http : / / presolar.wustl.edu/ pgd/). Model curves are from a 2 M© star of solar metallicity; 
only the final C-rich phases are plotted. Again, the model parameters are indicated. The 
range of carbon isotopic ratios of mainstream grains and of some A-l-B grains is reproduced, 
but our ^^N/^^N ratios are always much larger than solar. The full dots represent stellar 
measurements by Wannier et al. (1991) in C-rich circumstellar envelopes and confirm that 
high values are typical of evolved stars. See text for comments. The two symbols % and 
indicate the values of the nitrogen isotopic ratios found in Jupiter, which we use as initial 
abundance in our models, and in the Sun. 
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Fig. 13. — Left: a comparison of the NACRE rate for the reaction i'^N(p, 7)1^0 (hatched 
blue band) with the updated one from ADEIO (hght-red filled band). The solid lines mark 
the recommended value, together with lower and upper limits (blue for NACRE, red for 
ADEIO). The NACRE rate is adopted as a normalization (i.e. it is taken equal to 1 over 
the whole temperature region). Right: the same comparison for the reaction ^^N(p, 7)^®0. 
(A color version of this figure is available in the online journal.) 
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Fig. 17. — Same comparison as in Fig. 13, for the reactions ^^Mg(p, 7)^^A1, and 
25Mg(p,7)26Alg,. (A color version of this figure is available in the online journal.) 
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Fig. 18. — Same comparison as in Fig. 13, for the ^^Mg(p, 7)^^Al.m and ^^Mg(p, 7)^^A1 
reactions. (A color version of this figure is available in the online journal.) 
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